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We investigate the Rabi oscillation of an atom placed inside a quantum cavity where each mirror 
is formed by a chain of atoms trapped near a one-dimensional waveguide. This proposal was studied 
previously with the use of Markov approximation, where the delay due to the finite travel time of 
light between the two cavity mirrors is neglected. We show that Rabi oscillation analogous to that 
obtained with high-finesse classical cavities is achieved only when this travel time is much larger 
than the time scale that characterizes the superradiant response of the mirrors. Therefore, the 
delay must be taken into account and the dynamics of the problem is inherently non-Markovian. 
Parameters of interest such as the Rabi frequency and the cavity loss rate due to photon leakage 
through the mirrors are obtained. 


I. INTRODUCTION 


Many applications in quantum information process¬ 
ing rely on the strong interaction between stationary 
two-level emitters and photons [1]. This can be real¬ 
ized in cavity-QED (CQED) where an atom is placed 
inside a high finesse cavity [2-6], or waveguide-QED 
where the atom is placed inside or near a one-dimensional 
(ID) waveguide such as hollow-core fibers, fiber-taper 
waveguides and photonic crystal waveguides [7-12]. As 
CQED has become a well-established field, important 
experimental progress has been made toward enhanced 
coupling between atoms and the field in ID waveg¬ 
uides [ 11 , 12 ]. 

A hybrid strategy that combines the appealing features 
of both approaches has been proposed by Chang et al. 
in Ref. [13], where a chain of atoms trapped near a ID 
waveguide is used as a mirror in a cavity setup (see 
Eig. 1). In this seminal work, it is assumed that the col¬ 
lective response time of each atomic mirror is much longer 
than the time it takes for the photon to travel from one 
mirror to the other, and as a consequence the delay due to 
this travel time is neglected. This approach is the Markov 
approximation commonly used to accurately describe the 
interaction of photons with many atoms in a wide range 
of experimental situations. There are, however, certain 
cases where this assumption is not necessarily justified 
[14, 15]. Here, we show that in the Markovian regime the 
lifetime of the photon inside the cavity is not enhanced 
by the presence of the mirrors, and therefore sustained 
Rabi oscillation analogous to that observed in conven¬ 
tional CQED setups cannot be obtained. Indeed, it is in 
the non-Markovian regime, where the response time of 
the atomic mirrors is much smaller than the delay, that 
Rabi oscillation in the usual sense is achieved. 

As depicted in Eig. 1 , the problem under consideration 
involves an atom, initially in the excited state, located be¬ 
tween two quantum mirrors each consisting of N equally 
spaced atoms trapped near a ID waveguide. The posi¬ 
tion of this central atom is at an antinode of the stand¬ 


ing wave in the cavity in order to maximize the coupling. 
We study how the excitation probability amplitude of 
the central atom co(t) evolves with time. Our approach 
is valid for both Markovian and non-Markovian regimes. 
The figures of merit such as the Rabi frequency and the 
cavity loss rate through the mirrors are computed. Our 
main result is that the vacuum Rabi oscillation of the 
central atom is given by 


_ jt _ 

Co(t) e 2(l + iV7d/t;g)2 

where d is the distance between the two atomic mir¬ 
rors, 7 the single-atom decay rate into the waveguide 
modes and Vg the group velocity of light in the cavity. 
One sees from the above expression that the cavity loss 
rate is /^ = 7/(1 + N'^d/vgY and the Rabi frequency is 
^Rabi = 7y2V/(l +N^dlvg). 

The results obtained in Ref. [13], where flRabi ~ 
and the loss rate is comparable to 7, are retrieved in the 
limit N^d/vg 0, which corresponds to the Markov 
approximation. At the other extreme when d/vg ^ 
1 /(A/' 7 ), the lifetime of the photon inside the cavity is 
greatly enhanced similarly to conventional CQED se¬ 
tups. Since is the collective decay rate of the atomic 
mirror, 1/{N^) is the time scale of the mirror’s re¬ 
sponse. Thus, the strong interaction regime analogous to 
that achieved with high-finesse classical mirrors is inher¬ 
ently non-Markovian. In this regime, the Rabi frequency 
^Rabi ~ 7 ^ 2 i;^/( 7 (i) depends on the cavity length in¬ 
stead of the number of atoms in the mirrors. The ^1 /\fd 
dependence is expected because in one dimension d plays 
the role of the modal volume of the field in the cavity. 


2N 


1 + Njd/v^ 
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(1) 


II. QUANTUM CAVITY 

Let US first describe the atomic mirror formed by 
a chain of A" 1 identical two-level atoms strongly 
coupled to a one-dimensional waveguide and equally 
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FIG. 1. (color online). Cavity QED using atomic mirrors, 
(a) Scattering of an incoming photon on a chain of atoms with 
lattice constant dm- When uJAdmIvg — Itt with I an integer, 
the chain forms a Bragg mirror, (b) An initially excited atom 
(orange) is sitting inside a quantum cavity formed by two 
atomic mirrors (red). The spatial extension of the mirrors is 
small compared to the cavity length d. 



My 

FIG. 2. (color online). Reflectance of an atomic mirror 
as a function of the normalized detuning A/y. The green, 
blue and red line represent respectively a mirror made of 
N — 1,10, and 100 atoms. For reference, the dashed line 
illustrates the case where each atom forming the mirror fluc¬ 
tuates around its average position Xj with {dx‘^) = 0.01dm- 
iOAdmIc = TT. 


spaced by a distance dm^ as illustrated in Fig. la. For 
a monochromatic single photon impinging on the atomic 
mirror, the reflectance Rm strongly depends on the lat¬ 
tice constant dm and the detuning A between the fre¬ 
quency of the input photon and uja [16]. When the phase 
accumulated between two atoms of the chain ujAdm/'^g is 
a multiple of tt, the reflectance for a photon close to res¬ 
onance takes the form of a broadened Lorentzian [13] 

^ l + (A/(iV7))2- 


configuration the central atom is located at an antinode 
of this mode, maximizing its interaction with the cavity 
electromagnetic field. All the atoms are assumed to be 
strongly coupled to the ID waveguide and we first neglect 
losses due to the coupling with the environment. Effects 
of this noise will be discussed later. 

The atomic transition frequency uja between ground 
l^f) and excited |e) states is assumed to be much larger 
than the cutoff frequency of the waveguide. The dipole 
Hamiltonian describing the coupling between the atoms 
and the light field, under the rotating wave approxima¬ 
tion, is then given by [ 20 ] 


The reflected frequency bandwidth is thus which is 
significantly enhanced by the large number of atoms in 
the chain, as shown in Fig. 2. Moreover, while the re¬ 
flectance is sensitive to small fluctuations in the atomic 
positions Xj + 5xj for some specific frequencies, the over¬ 
all reflection coefficient of a broad bandwidth photon re¬ 
mains close to unity. This is an example of the well- 
known Bragg mirror [17, 18] and we will focus on this 
geometry for the rest of the paper. The presence of N 
atoms in the mirror also leads to a superradiant decay at 
a much higher rate A/'y, as recently observed experimen¬ 
tally [19]. The time scale of the collective response of the 
atomic mirror is thus which will play a key role 

in the following. 

We now consider the complete quantum cavity system 
formed by a central atom located between a set of two 
atomic mirrors [13], as illustrated in Fig. lb. This central 
atom is identical to those of the mirrors. An important 
parameter controlling the interaction between the central 
atom and the cavity is the phase shift acquired by a reso¬ 
nant photon across the cavity 0 = ujAdIvg. We will focus 
on the case 0 = ( 2 n+l) 7 r, where n is an integer, such that 
uja corresponds to a natural mode of the cavity. In this 


'dip = -ih ^ [ dio [4 

n —_ AT ’^0 

. ( 3 ) 


j=-N^ 


-H.c. 


where g^; is the coupling constant, Xj is the position of 
atom j with raising ladder operator = \e)j{g\, and 
and buj respectively annihilate right-going and left-going 
photons. In this work we use the Weisskopf-Wigner ap¬ 
proximation and set g^; = goo a- Under this assumption, 
the decay rate of a single atom into each waveguide spa¬ 
tial mode is 7 = 27r^^^ [20]. It should be noticed that 
under typical experimental conditions, the cavity oper¬ 
ates in a regime where ^djvg 1 N (see Table I). 
In particular, in the limit N ^ oo where the mirrors 
would ideally reflect all frequencies, one would expect to 
recover Rabi oscillation between the central atom and 
the cavity field as described in the framework of cav¬ 
ity quantum electrodynamics. However, as emphasized 
in the introduction, the key parameter which determines 
the dynamics of the quantum cavity is given by the ratio 
between half the cavity round-trip time and the atomic 
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TABLE L Experimental parameters for a trapped cesium 
atom [19] and quantum dot [11] coupled to a photonic-crystal 
waveguide, and superconducting qubit coupled to a ID copla- 
nar waveguide transmission line [21]. Eor the distance be¬ 
tween the mirrors we consider the range of values usually en¬ 
countered in conventional CQED experiments with each type 
of TLEs. 


TLEs 

a;A(GHz) 

27(MHz) 

27 / 70 "" 

Vg/C 

d(mm) 

'ydjvg 

Cs atom 

2.1 X 10“ 

32 

1.1 

0.1 

1 

5.3 X lO"'^ 

QD 

2 X 10® 

6.2 X 10® 

63 

0.01 

10“^ 

1.0 X 10"^ 

SC 

7.1 

6 X 10^ 

> 20 

0.5 

10 

2 X 10"^ 


^ 70 is the decay rate into the environment, i.e. outside of the 
waveguide modes. 


mirrors response time N^d/vg. 


III. TIME EVOLUTION 


We consider the situation when the atoms in the mir¬ 
rors are initially in the ground state, while the central 
atom is prepared in the excited state. The Schrodinger 
equation yields the following delay differential equations 
for the atomic excitation amplitudes (see Appendix A) 


co{t) =-'yco{t) - 
Cmit) = - 7V^e®^/^co(i - 


jN Cm{t) + Cm{t - ^)Q{t - ^) 


(4) 

where we have formally integrated the field variables, 
~ yiw depends on the excitation 

amplitude of the atoms forming the mirrors and 0(t) is 
the Heaviside-step function. Here I is the integer defined 
by ujAdm/'^g = Moreover, the delay due the traveling 
time of a photon exchanged between the central atom 
and the two mirrors forming the cavity is included. 

These equations can be solved by Laplace transform. 
It is straightforward to show that the transform of co(t) 
is 


smaller than the collective response time 1 / {N^) of the 
atomic mirrors. By setting d = 0 in Eq. (5) and obtain¬ 
ing the inverse Laplace transform, we find that, when 
V » 1, co(t) = cos(V 2 V 7 t), which is consis¬ 

tent with the results of Ref. [13]. In other words, the 
central atom under goes damped oscillations at the fre¬ 
quency r^Rabi = V 2 V 7 , which however is not similar 
to the sustained Rabi oscillations often encountered in 
conventional CQED. Indeed, the associated decoherence 
rate 7/2 is found to be comparable to the decay rate of a 
single atom into the waveguide modes when the mirrors 
are not present, which suggests that the photon is actu¬ 
ally not trapped in the cavity. Moreover, the oscillation 
frequency depends on the number of atoms in the 

mirrors, while the Rabi frequency usually depends on the 
modal volume of the field inside the cavity [ 22 ] , which for 
our ID architecture is given by the distance d between 
the two mirrors. 

To provide further insight on why the oscillations in 
the Markovian regime do not show the usual features 
observed in conventional CQED, we recall that by ne¬ 
glecting the delay one has assumed that the central atom 
and the mirrors feel the influence of each other instan¬ 
taneously, i.e. d/vg <C 1 /( 7 ^). Now for a propagating 
photon to be trapped in the cavity, its real-space distri¬ 
bution [23] must have a width at most comparable to d. 
This implies that the spread of this photon in momentum 
space is bounded by Auj > Vg/d. On the other hand, Auj 
has to be smaller than the reflection bandwidth 7 V of 
the atomic mirrors, which implies that d/vg > 1 /( 7 ^), 
where 1 / (jN) is the collective response time of the mir¬ 
rors. This inequality manifestly contradicts the Markov 
approximation which assumes the delay to be smaller 
than any other dynamical time scale in the cavity. 

The above discussion makes it clear that the dynamics 
of the system is inherently non-Markovian in the regime 
where the photon is trapped in the cavity. Moreover, the 
critical parameter that determines the transition between 
the two regimes is N^d/vg. 


B. Non-Markovian regime 


g ^ 5 + 7A^(1 + 

(s + 7)(s + 7A/') + — 7)7A/' 

The excitation amplitude cq (t) is then obtained by taking 
the inverse Laplace transform. We show below that a 
compact analytical expression for co{t) is possible in the 
regime of interest jd/vg 1 N. As mentioned above 
we will mainly consider the case 0 = (2n + l) 7 r, which is 
assumed in the rest of the text unless stated otherwise. 


A. Markovian regime 

We first apply the Markov approximation and ne¬ 
glect the delay d/vg, which is valid if this delay is much 


We now study the time evolution of the central atom in 
the non-Markovian regime, where N^d/vg may not nec¬ 
essarily be small. This requires solving Eq. (4) without 
neglecting the delay. The typical procedure for obtain¬ 
ing the inverse Laplace transform is to expand co(s) in 
a geometric series and then taking the inverse transform 
of each term in the series [24]. As shown in Appendix B, 
this yields 

00 

co(i) = kd/vg)Q{t - kd/vg), (6) 

fe =0 

where fk{t) are smooth functions in time, /o( 0 ) = 1 and 
//c>i(0) = 0- The exact form of fk{t) is not needed for 
our discussion and can be found in Appendix B. The 
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relevant observation here is that the Heaviside function 
gives rise to discontinuity in the slope of co(t) at times 
kd/vg^ which correspond to the retardation times asso¬ 
ciated with successive exchanges of photons between the 
central atom and the mirrors (see Fig. 3d). This behavior 
is often encountered in the dynamic of quantum systems 
with delay [25, 26]. 

However, under typical experimental circumstances, 
^djvg C 1 , and thus the evolution of the central atom 
appears smooth when the time scale of observation is 
comparable to I/ 7 . With this condition, it is also pos¬ 
sible to obtain a compact analytical approximation for 
Co (t) which gives much more insight into the dynamics of 
the central atom than the exact series solution in Eq. ( 6 ) 
does. In order to see this more clearly, we first discuss the 
“macroscopic” regime when Njd/vg ^ 1 so that we can 
set ^ 00 in co(s)^. As a result, the Laplace transform 
in Eq. (5) is simplified to 


Co{s) 



1 + e 


1 — e 


SdjVg 

SdjVg 


-1 


( 7 ) 


Instead of geometric series expansion, another ap¬ 
proach for obtaining the inverse Laplace transform is 
to identify the poles Sj of co(s). The excitation ampli¬ 
tude is then readily obtained via Cauchy’s theorem as 
co(t) = Sj where Vj is the residue of co{s)e^^ evalu¬ 
ated at the pole Sj. We find that the poles of co(s) are 
simple poles confined to the imaginary axis, i.e. sj = iyj 
where yj is a real number. Moreover, the exact location 
of the poles are given by the solutions of the equation 
y = 7 cot(^), which leads to an infinitude of poles that 
have reflection symmetry through the coordinate origin. 
The residue at the pole sj is rj ~ e^7^/[l — s^d/{2^Vg)]. 
When jd/vg 1, the two poles nearest to the origin 
are s±i ^ ±iy^ 27 c^/d, and the contribution of these two 

poles to co(t) is approximately cos (^^2^Vgld The 
contribution of the other poles are shown to be bounded 
above by ^d/{6vg), which is negligible (see Appendix B). 
Therefore, we obtain for the central atom excitation am¬ 
plitude 


Co{t) ~ cos 



(8) 


This results implies that in the macroscopic regime 
N^d/vg ^ 1 , the central atom undergoes sustained Rabi 
oscillation without decoherence. This is understandable 
since the atomic mirrors have a reflection bandwidth that 
is much larger than the frequency width of the photon 


^ Although N ^ Vg/{'yd) in this regime, the assumption that the 
spatial extension of the mirror is small compared with the dis¬ 
tance between the mirrors is still valid as long as Ndm d. We 
can thus safely neglect any increase in the cavity mode volume 
due to penetration of the field into the mirrors [27]. 



FIG. 3. (color online). Probability of excitation of the cen¬ 
tral atom for different values of the parameter Njd/vg. The 
green curve illustrates the analytical approximation (10) de¬ 
rived in the regime 'ydjvg <C 1 ^ N, while the red dashed 
curve is obtained for any regime by directly inverting the 
Laplace transform in the form of (6). The decaying enve- 
__ 

lope e jg plotted in dashed black for reference, 

(a) In the Markovian regime, we retrieve the damped oscil¬ 
lation obtained in [13], where the photon is not trapped in 
the cavity, jdjvg = 0.0002 (b) In the transition regime, the 
lifetime of the photon in the cavity is significantly enhanced. 
'ydjvg = 0.01 (c) In the macroscopic regime, the photon is ef¬ 
fectively trapped in the cavity and we observe sustained Rabi 
oscillation between the central atom and the cavity mode field. 
'ydjvg = 0.02 (d) For a large cavity length ^djvg — 0.5, the 
retardation effects give rise to discontinuity at times yt = 0.5/c 
with k an integer. This is not captured by Eq. (10) which is 
limited to the regime 'ydjvg 1. 


trapped in the cavity. Moreover, the Rabi oscillation is 
proportional to 1 /\/d, which is expected since in one di¬ 
mension d plays the role of the modal volume of the field 
in the cavity. 

Now we consider the more general situation when 
Njd/vg may have arbitrary value. The two poles nearest 
to the coordinate origin are shifted to 


S ±1 « -7/(2(1 + N'jd/Vg)^) ± + Njcl/Vg), 

( 9 ) 


and their contribution yields 


co(i) 


_ Jt _ 

g 2(l + N^d/vg)^ COS 


(\/i+7/../)' 


( 10 ) 


The contribution of the other poles are again negligi¬ 
ble in the regime of interest ^djvg ^ 1 ^ N. Both 
the Markovian limit, Njd/vg 1, and the macroscopic 
limit, Njd/vg ^ 1, are correctly retrieved from Eq. (10). 
Eig. 3 shows a comparison between our compact analyt¬ 
ical approximation obtained in Eq. (10) and the exact 
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form of co{t) in Eq. ( 6 ) for various values of jdjvg and 
N. The approximation is indeed very accurate in the 
regime of small ^djvg. 

From Eq. (10), one sees that the cavity loss rate is given 
by/^ = 7 /(l + N^d/vgY while the Rabi frequency reads 

^Rabi = 1\J define the critical size of the 

mirror Nc = Vg/{jd) that signals the transition from 
the delocalization to the localization of the photon in the 
cavity. The cavity loss rate is greatly suppressed when 
N > Nc. One can refer from Table I that Nc ~ 1900 for 
the first configuration with Cesium atoms, Nc ~ 100 for 
quantum dots and Nc ~ 50 for superconducting qubits. 
While the first configuration seems to require a large 
number of atoms, it should be noted that the critical size 
Nc can be reduced by an order of magnitude by slow¬ 
ing light in the photonic-crystal waveguide [28, 29] or 
increasing the coupling strength between the atoms and 
the waveguide. Besides, as discussed previously, small 
position fluctuations do not change much the average re¬ 
flection of the mirrors forming the cavity; therefore, one 
does not need a significant extra number of atoms in or¬ 
der to trap the photon in the cavity. 

IV. EFFECTS OF DETUNING AND LOSS TO 
ENVIRONMENT 

We have studied the Rabi oscillation of the central 
atom inside a quantum cavity when the atom is on reso¬ 
nance with one of the natural modes of the cavity, that 
is, uja = (2n + X^'KVgjd. We now study the effect of a 
small detuning between the atomic transition frequency 
and the cavity modes and show that our model recov¬ 
ers the behaviour predicted by CQED with classical mir¬ 
rors. In particular, the phase accumulated across the 
cavity by a photon on resonance with the atom now 
reads d = (2n + l)7r-h0 with (f ir. The frequency 
of the main cavity mode interacting with the trapped 
atom is cjc = (2n + l)7TVgld^ yielding a detuning A = 
uja — oJc = (fvg/d. Focusing on the macroscopic regime 
(see Appendix C for the case of arbitrary N^d/vg), the 
poles of co(s) in Eq. (5) are now given by the solutions 
of the equation y = 7 cot(^^^^^^). The two main poles 

are thus shifted to = iA /2 ± + (A/ 2 )^, where 

Uq is the Rabi frequency obtained previously without de¬ 
tuning, and their contribution yields 

co(t) ^ cos(Ut) — i sin(Ut)^^ ( 11 ) 

where the generalized Rabi frequency reads Q = 
a/Uq + (A/ 2 )^, similarly to conventional CQED setups 
[30]. 

In realistic circumstances, each atom coupled to the 
ID waveguide also emits independently into the free- 
space environment with rate 70 . In practice, this in¬ 
duces additional loss in the cavity and impacts the deco¬ 
herence rate of the Rabi oscillation obtained in Eq. (10) 


which would now read b: + 70 . We thus see that the 
decoherence rate is now bounded from below by the 
loss to the environment at 70 , no matter how the key 
parameter N^d/vg is optimized. However, waveguide- 
QED setups allow the realization of remarkably high 
coupling efficiency between the atoms and the waveg¬ 
uide modes, which implies that the period of Rabi os¬ 
cillation is much smaller than the time scale of the loss 
for large N. Indeed, for the three configurations in Ta¬ 
ble I, the ratio between ^Rabi and the total decoherence 
rate is around 21, 35 and 20 at the critical number of 
atoms Nc in the mirrors . Hence, one expects to observe 
many cycles of Rabi oscillation before the photon even¬ 
tually leaks to the environment. The cooperativity is 
given hy 7] = fJRabi/( 7 oK) = 2Ar(l + Nyd/c)y/yo, which 
is very large when N is close to Nc. Large cooperativ¬ 
ity is essential for many quantum information processing 
applications based on strong atom-photon coupling such 
as the nanophotonic quantum phase switch [31] and the 
atom-photon quantum logic gate [32]. 


V. CONCLUSION 

When comparing this work with the description of a 
cavity made of classical mirrors, the reader might be sur¬ 
prised by the apparent need of a non-Markovian theory in 
order to reproduce standard Rabi oscillation. However, 
one should notice that the non-Markovianity is required 
to build up the cavity standing wave out of propagating 
modes. This is not in contradiction with standard CQED 
where the standing wave structure is assumed a priori by 
imposing that the electric field vanishes at the mirrors po¬ 
sition. In the latter scenario, one then assumes that the 
evolution occurs simultaneously throughout the cavity, 
with the photons being delocalized in the standing-wave 
spatial mode. Notably, this description is valid in the 
regime ^djvg 1 where retardation effects are negli¬ 
gible [33]. A non-Markovian treatment is required only 
when we also consider the dynamics of the mirrors whose 
characteristic time scale is much smaller than the delay, 
as can be seen in our calculation. 

In this work, we have derived the Rabi frequency and 
the decoherence rate of the Rabi oscillation undergone 
by an atom placed inside a quantum cavity formed by 
atomic mirrors. We found a simple condition for testing 
whether the dynamics of the system can be correctly de¬ 
scribed with the Markov approximation. Our approach is 
valid for both the Markovian regime, where the photon’s 
travel time between the mirrors is neglected, and the non- 
Markovian regime, where we found that the photon can 
be trapped inside the cavity. In the latter regime, we 
have found that sustained Rabi oscillation, as achieved 
in conventional CQED with high finesse mirrors, can be 
observed and the experimental parameters required for 
achieving it are estimated. 
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Appendix A: Deriving the coupled delay-differential equations 


As illustrated in Fig. lb, our system is constituted of two chains of N atoms separated by a distance d, with a 
spacing of dm between neighboring atoms in each chain, and an additional central atom located between both chains. 
We consider the situation when the length of the chain is much smaller than the distance between the two atomic 
mirrors. The Hamiltonian in the interaction picture reads 


Hdip = Y] / - H.c. 

j=-JV -^0 

where is the position of the jth atom. The state of the system can be written as 

^ nOO 

I^W) = Y CjW'5'+|0)+ / dLo{ca{co,t)al + Cb{Lo,t)bl)\0), 

j=-N “'o 


(Al) 


(A2) 


where |0) indicates the combined atom-field vacuum, which is the state when there is no photon in the waveguide and 
all the atoms are in the ground state. The initial state where only the central atom is in the excited state corresponds 
to Cj{0) = Sqj and Ca(ci;,0) = C5(co’,0) = 0. The Schrodinger equation then leads to the following set of equations 


N 


Ca{(x,t) = Yi 

-iOJXj /Vg 

5 

1 

II 


N 


Cb{u;,t)= Y 

j=-N 


nOO 

Jo 



(A3) 

(A4) 

(A5) 


In the Weisskopf-Wigner approximation one can replace by which is done in the rest of the text. Integrating 
formally the differential equations for Ca{uj,t) and and then inserting them into the equations for Cj(t), we 

obtain the following delay-differential equations for the atomic excitation amplitudes 


N 

Ci W = -7 E - \x^ - x,,\/v,). 

j' = -N 


(A6) 


For the configuration of our system, the atoms are located at the following positions 


{ xo = 0, 

Xj = {j -h l)dm - d/2 for - N <j < -I, 

Xj = {j - l)dm + d/2 for I < j < N. 

Since the spatial extension of each atomic mirror is much smaller than the distance between the two mirrors, we 
can neglect the time delay due to the distance between any two atoms in the same mirror. Moreover, we focus on 
the configurations where ooAdm/xg = lir for two adjacent atoms in the same mirror, where I is an integer. In this 
scenario one can check that every coefficient { — for j ^ 0 are equal. Denote the excitation amplitude of 
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the cavity by Cm{t) = we then obtain the two coupled delay-differential equations for the 

central atom and the cavity 


Co{t) = -JCo{t) - jVm- ^)0{t - 2^), 

Cm{t) = - 2^)e(t - 2^) - 7iV [emit) + - ^)0(i - 

For 9 = (2n + l)7r, the Laplace transform of co(f), denoted by co(s), is 

s + 7 A ^(1 - 1 

^ “ is + 7 ) (s + JN) 1 _ e-^d/v, 1^7)7^ ’ 

(s+7)(s+7Ar) 

and Co (t) can be obtained by inverting the Laplace transform, which in our case can be done by summing the residues 
at the poles of e^^co(s). 


(A 7 ) 

(A8) 

(A 9 ) 


Appendix B: Solving the equation 

Let us expand the second term in the expression of co(s) in a geometric series 


1 


I _ ^-sdlvg is-^)-fN 


= E 


“ " (s-7)7Ar 


^—skdlvg 


(s+7)(s+7A/') k=0 


is + 7)(s + jN) 


We then get 


Co(s) = 


s + 7 A 


(s + 7)(s + 7Af) ^ 


E< 


—skd/vg 


-fN 


{s -h 7)(5 + 7A) 


k=0 


E« 

k=0 


{s — j)jN 


1 k 


{s + 7)(5 + 7 A) 


— s{k-\-l)d/vg 


{s — 7)7^ 

{s + 7)(s + 7 A) 


n k 


and by taking the inverse Laplace transform term by term we have 

00 

^od) = '^fkit- kd/vg)Qit - kd/vg) 


/c=0 


(Bl) 


with 


= e 


fo = ^-"{^} = 

I 5 + 7 J 

fk>i = 


-7t 


r is - 7)777 

* 777 

(s — j)jN 

k-1 ^ 

[(s + 7 )(s + 777). 

(s + 7)(s + jN) 

_(s + 7)(s + 'yN)_ 

I 


= c 


-'{I 


{s — j)jN 


1 k 


+ 7 s-'jl [(s + 7)(5 + 7A)_ 


where C~^ denotes the inverse Laplace transform and we made use of = C~^{F}{t — kd/vg)Q{t — 

kdjvg). One can check that fk is a smooth function of time since its Laplace transform involves only polynomials, and 
furthermore //c( 0 ) = 0 for /c > 1 . The evolution for t < d/vg is that of a single atom in the waveguide, which is due to 
causality: For t < dji^Vg)^ the photon wavepacket emitted by the central atom has not yet reached the mirrors; for 
t > dl{2vg\ the mirrors can be excited by the inffuence of the central atom, however their back-action on the central 
atom will be delayed by another duration of dj(2vg). 

The series solution for co{t) given in (Bl) is exact but does not give much insight into the dynamics of the central 
atom. We now derive the compact analytical approximation of co(t) in the regime of interest jdjvg <^1 N. First, 
in the Markovian regime when Njd/vg 1 , the solution is straightforwardly obtained by setting d/vg = 0 in Eq. (A 9 ) 
and taking the inverse transform, which yields 


Cod) 


_ ji 

= e 2 cos 




O 



(B 2 ) 

































)/ 6 l{ 2 Vg) 


FIG. 4. (color online), y' in red and coi{y')'yd/{2vg) in blue as a function of the dimensionless variable y' = ydl{2vg), for 
^djvg = 0.1. The intersections correspond to the poles of co(s). 


In the opposite extreme, the “macroscopic" regime when Njd/vg ^ 1, we set ^ oo in Eq. (A9) and obtain 


Co{s) 


I _ ^-sdjvg 

s + 7 — (s — 7)e“^^Ag 


(B3) 


The poles, which are given by the zeros of the denominator, must obey = e ^ and thus = e ^ 

With a bit of algebra one can show that this constraint is satisfied if and only if Re(5) = 0. Hence the poles of co(s) 
must lie on the imaginary axis. Denote y = Im(5), we have 


CQ{iy) = 


—i 

2/-7cot(|J)’ 


(B4) 


and the poles are given by the solution of ^ = 7 cot(^). Note that both functions are odd, so if ^ is a solution 
then —y is also a solution (see figure 4). In the limit ^djvg ^ 1, the two solutions closest to the origin are given by 

y±\ There are also other solutions close to each singularity of cot(^), denoted by = ±2, ±3,- 

From Fig. 4 we see that \yj\ > {\j\ — l)27TVgld. 

We can now obtain co(t) by evaluating the residue at each pole 


co{t) 


The sum can be bounded as follows 


w lim — y-y^ ^ , 

^V^Vj y-jcot{§^) 


dyt 


E 

j 




1 + J± + J^Lcot('^'l2 

^ ^ 2v^ ^ 2v^ cuuu 




2Vg 

Ay-it 


1 I ^ 


2 d 
y±l 2-iVg 


y2vg 

oo 


^ 1 _L 

i=2 ^ 2v, 


givjt _(_ giy-jt 


■ ^ 
^3 ‘^IVg 


E 


i=2 


2 cos{yjt) 


1 I ^ 


y 


2 d 

'j 2-iVg 


oo 


<E 


^3 2^Vg 


1 oo 

<^y _ 


1 yd 

6Vg 


(B5) 


(B6) 


since \yj\ > {\j\ — V)2'KVgld for j ^ ±1. Therefore, we arrive at the analytical approximation 
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One observes that the main contribution to co(t) is from the two poles nearest to the coordinate origin. We find 
that there is a simpler procedure to obtain the position of these two poles approximately. This is done by expanding 
e~sdlvg ^ l—sdlvg-\-{sdlvg)‘^l2 in the expression of co(s) in Eq. (B3) and finding the zeros of the resulting denominator 
(which is a cubic polynomial but only two of its zeros correspond to the two main poles we are interested in). 

Now we study the general case when N^djvg = a where a can have any finite value. First replace N = avg/{dj) 
in the expression of co{s) in Eq. (A9), and then expand to second order. Solving for the zeros of the resulting 

denominator, we find that the two main poles are now shifted to 


5±i 


-7/[2(1 + a)2] ±i7 


2avgl{d‘y) 
1 T (2 


and their contribution yields 


_ 7* 

Co(t) = e 2 (l + a)2 






(B 8 ) 


(B9) 


which is Eq. (10). 


Appendix C: Effects of Detuning 


Let us consider a detuning A between the cavity and the atomic frequency. Then 0 = (2n + l) 7 r + 0, where 
(j) = Adjvg. The Laplace transform now takes the form 


Co{s) = 


s + 7^(1 - 


{s + 7)(5 + 7 A) 


1 piAdjvg f^ — sdjvg 7 ) 7 ^ 

^ ^ ^ (s-h 7 )(s-hjN) 


(Cl) 


We first focus on the good cavity limit when the photon is trapped, which happens when Njd/vg ^ 1 . After 
setting N ^ 00 we have 


Co{s) 


2 _ ^iAdjvg^-sdjvg 
S 'J — (s — ry'^^iAdjvg ^—sdjvg 


(C 2 ) 


As in the case without detuning, the poles lay on the imaginary axis and satisfy y = 7 cot(^^ 2 ^^^^). Hence, for 
Adjvg IT, the two main poles are now given by 


s±i = iA/2 ± + (A/2)2 


(C3) 


with flo = Repeating the derivation of (B5), the contribution of the other poles is still of order 0{^^djvg) 

and the inverse Laplace transform reads 


co(t) = e* 2 ^ ^ cos(flt) — i sin(Ot)^^ + 0{'ydfvg) 

which is consistent with cavity-QED. The probability of excitation is then 

cos(m) + (^)^ 


Po{t) = 




(C4) 


(C5) 


For the general case of arbitrary detuning and arbitrary value of the parameter a = Njd/vg, we also found an 
analytical approximation that matches the exact formula in Eq. (Bl) to a high accuracy when jd/vg 1 which 

is given by 


co(t) = exp ( - 


7t 


2(a + l)2 


exp 


2a{a + l)v 
SjdjVg 


7t 


cos (floVl -^u^vt) 


+ 


1 ^ U ^ V 


sin (floVl -^u^vt) 


0{jd/Vg) 


(C 6 ) 
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FIG. 5. (color online). Evolution of the excitation of the central atom, for different regimes, with ^djvg ^ 1 ^ N and 
0 = tt/IO. The exact solution from Eq. (Bl) is shown by the solid blue line, and the analytical approximation from Eq. (C6) 
the dashed red line, a) 'ydjvg = 0.0002. b) ^djvg — 0.01. c) 'yd/vg — 0.02. 


with llo = ~ ^ ~ f J. The strong agreement between this formula 

and the exact series solution for different values of a corresponding to the Markovian, non-Markovian and transition 
regimes can be seen in Fig. 5. 
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